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Abstract. We employ the very recent photometric data 
of the OGLE project together with stellar atmosphere and 
linear pulsation models to determine the distance modu- 
lus of the Small Magellanic Cloud from its double-mode 
Cepheids. Based on the requirement of obtaining the same 
distance modulus DM from the two types of variables 
(fundamental & first overtone and first & second over- 
tone), we get DM = 19.05 mag, with a very small sta- 
tistical error (standard deviation) of 0.017 mag. Various 
systematic and zero point ambiguities (primarily those of 
the color-temperature transformation) lead to an error of 
±0.13 mag (estimated 3cr deviation). This result is in very 
good agreement with the distance modulus of the Large 
Magellanic Cloud of 18.5 mag, derived earlier from cluster 
double-mode RR Lyrae stars. 
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1. Introduction 

In a former paper (Kovacs & Walker 1999, hereafter 
KW99) we have shown that the distances derived from 
double-mode RR Lyrae (RRd) stars are systematically 
larger than the ones obtained from the standard Baade- 
Wesselink (BW) analyses of RR Lyrae stars. The latter 
distances are in similar conflict also with the BW results 
of Cepheids (e.g., Gieren et al. 1998). The reason of this 
discrepancy is still unknown (Cacciari et al. 2000). This 
contributes further to the ambiguity at the 0.2-0.3 mag 
level over the luminosity scale of the RR Lyrae stars, and 
consequently, over the distances to the nearest globular 
clusters and galaxies. 

The purpose of this Letter is to study the applicabil- 
ity of double-mode variables in the distance calibration in 
more detail. The discovery of a large sample of fundamen- 
tal & first overtone (FU/FO) and first & second overtone 
(FO/SO) Cepheids in the Small Magellanic Cloud (SMC) 



by the OGLE team (Udalski et al. 1999a, hereafter U99) en- 
ables us to derive the distance to the Cloud and thereby 
adding another piece of information to the dispute over 
the distance of the Magellanic Clouds. 

2. Data, method, models and temperature scales 

Double-mode (or beat) Cepheids in the SMC have been 
discovered previously by the major microlensing projects 
(macho, Alcock et al. 1997; eros, Beaulieu et al. 1997). 
However, it is only the OGLE team who publishes the data 
in standard colors (i.e., in Johnson V and in Kron-Cousins 
I c ). Therefore, for the time being, we decided to employ 
their data only. 

We use the periods, average V — I c colors and V magni- 
tudes of the 23 FU/FO and 70 FO/SO Cepheids published 
by U99{|] Although they presented also B — V colors, we 
decided not to use them, because of the few data points 
in B. For reddening we accepted their values derived for 
the various fields from a method based on red clump stars 
(see U99 and the ftp site mentioned above) . 

In computing the distance modulus, we follow almost 
entirely the method of KW99. Here we repeat only the 
basic steps and assumptions. 

For any double-mode variable, from the pulsation mod- 
els we obtain relations between the physical parameters 
and the periods 



Pi = fi{M,L,T^,X, Z) 



(1) 



where i = 0, 1 or 1, 2 for the FU/FO and FO/SO variables, 
respectively. The other parameters have their usual mean- 
ing. In principle, from the observed color we can determine 
T c g and from other information we also have approximate 
values for the hydrogen and metal abundances. Therefore, 
we can invert the above relations to compute the luminos- 
ity (and mass). Next, the distance modulus DM is com- 
puted simply from the comparison of the observed V mag- 
nitude and the calculated absolute magnitude from L by 
using a bolometric correction (B.C.) formula and proper 
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interstellar reddening. The existence of the two types of 
beat Cepheids in SMC enables us to estimate DM in two 
different ways. Because of the different dependence of the 
periods on the physical parameters in the two types of 
variables, this will allow us to calculate DM in an op- 
timum way and put constraints on Z, one of the most 
important parameters entering in this method. The other 
crucial parameter is the zero point of T e g, which cannot 
be constrained from the present data, but, in principle, it 
is also possible to do. 

The main assumptions entering in our approach are 
the following: 

— FU/FO and FO/SO variables have the same chemical 
composition. 

— Linear nonadiabatic and purely radiative pulsation 
model periods are applicable to the observed periods. 

— The SMC has small spatial extent relative to its dis- 
tance. 

— Current color-temperature calibrations are reliable 
enough. 

As regards the technical details, it should be mentioned 
that we used the same standard, purely radiative linear 
pulsation code as in KW99. A large number of models 
were computed with X = 0.76 and Z = 0.001, 0.002, 
0.003, 0.004, 0.008 and using OPAL'96 opacities (Iglesias 
& Rogers 1996). The distribution of the various elements 
in Z corresponds to that of the Sun. Unlike in the case 
of RRd models, we could not find a simple way to invert 
Eq. (1). Therefore, in finding the values of (M,L) which 
fit the observed periods the best for the given T e g, X and 
Z, we used a straightforward search in the interpolated 
fine grid of the original model sequences. 

As far as the (V — I c ) — ► T e g and B.C. calibrations 
are concerned, we used the stellar atmosphere models of 
Castelli et al. (1997) and obtained the following formulae 
through least squares fits in a parameter space relevant for 
beat Cepheids (i.e., T eff = 5000-7000K, log g = 2.0-3.5, 
[M/H] = -1.5-0.0) 

logoff = 3.9224 - 0.2470(U-/ C )+ 0.0046 log 5 

+ 0.0012[M/H] , (2) 
B.C. = 0.0411 + 2.0727AT- 0.0274 \ogg 

+ 0.0482[M/H] - 8.0634AT 2 , (3) 

where AT = logT e ff — 3.7720, and we used the stan- 
dard notation for the relative heavy element abundance 
[M/H]= logZ/Z Q , with Z Q = 0.02. 

The above B.C. is adjusted to M bo i(0) = 4.75 and 
yields B.C.(q) = -0.09. The zero point of Eq. (2) is 
only marginally higher than that of Blackwell & Lynas- 
Gray (1994, hereafter BLG94) (see also Clementini ct 
al. 1995), which is based on the InfraRed Flux Method 
(IRFM). Other, more recent IRFM-based T off calibra- 
tions (Blackwell & Lynas-Gray 1998, hereafter BLG98; 



Alonso et al. 1999, hereafter A99) yield lower zero 
points (log T eff (BLG9&) - logT off (BLG94) w -0.004, 
logT eff (A99) - log T eff (5LG94) w -0.010). We caution 
however, that these differences are based on (V — I c ) — ► 
T e g calibrations, and the corresponding formula of A99 
can be employed only after applying a transformation 
from their Johnson I to I c (Fernie 1983). This could per- 
haps be one of the reasons why we get (somewhat curi- 
ously) better agreement between BLG94 and A99 with the 
temperatures calibrated by B — V. In the next section we 
will check the sensitivity of the derived distance modulus 
against the various T e g scales. 

Finally we note that for the calculation of the gravity 
we used the following formula 

logg = 2.62- 1.21 log P FU , (4) 

which has been derived from a simple pulsation equation 
and black-body relation. This formula has a maximum 
error of ±0.05, assuming that M/M Q = 3.0t?;g. 

3. Determination of the distance modulus 

To characterize the quality of the fit, for each variable, we 
calculated the following quantity 

a = ^Aif + A/f +1 , (5) 

where i = or 1 for the FU/FO and FO/SO variables, 
respectively. The difference between the observed and cal- 
culated periods are denoted by APj. We found basic dif- 
ference between the two types of variables in respect of 
the behavior of a(M,L). In Figs. 1 and 2 we show rep- 
resentative gray maps and one dimensional slices along 
the minimum values of a. For better visibility, in the gray 
maps we used a relatively low resolution, and therefore we 
considered only the minimum value of a in each pixel (the 
scans were performed on a much finer grid). We see that 
the minimum (indicated by the black pixels) is much more 
pronounced for the FU/FO than for the FO/SO variables. 
This is nicely seen also in the slices along the minimum 
values. In addition, there is also a slight offset of < 0.002 
for PijP\ between the observed and theoretical values. 
Although this problem might bear some theoretical sig- 
nificance, it is completely unimportant in the present con- 
text (see Table 1). Furthermore, the shallow minima of 
the FO/SO variables yielded more stable DMs, depend- 
ing (much) less on the various parameter changes than 
those of the FU/FO variables. This, together with their 
large number, make them very valuable for the purpose of 
distance estimation. 

The calculation of the distance modulus was performed 
for each Z listed in the previous section. By applying the 
temperature scale given by Eq. (2), the individual DMs 
are plotted in Figs. 3 and 4. We employed the 3a crite- 
rion for filtering out a few outliers (the number of these 
variables never exceeded two). The most striking feature 
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Fig. 1. Lower panel: Gray map of the period deviation 
(Eq. (5)) in the (M, L) parameter space. Darker areas in- 
dicate better fits to the observed periods. Upper panels: 
minimum period deviations as functions of M and log L. 
An FU/FO variable is tested. 
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Fig. 3. Individual and average distance moduli calculated 
for the FU/FO and FO/SO variables (shown by open and 
gray filled circles, respectively). Average DMs are shown 
by thick lines. Chemical compositions of the models used 
are given in the top left corners. 
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Fig. 2. As in Fig. 1, but for an FO/SO variable. 



of these plots is the large sensitivity of DMq\ (DM of 
the FU/FO variables) against the variation of Z. This 
is in contrast with the very weak sensitivity of DM\2- 
These properties and the DMqi = DM\2 condition al- 
low us an optimum estimation of Z . (We note in passing 
that the other natural criterion, minimum dispersion of 
all individual distance moduli, leads to the same conclu- 
sion.) As it is seen in Fig. 4, the present results prefer a 



somewhat lower abundance than the one most often used 
in the context of the SMC Cepheids. With Z Q = 0.02, 
our close to optimum Z corresponds to [Fe/H]= —0.8, 
whereas the usually quoted value is —0.7 (corresponding 
to Z = 0.004, see Luck et al. 1998). We think that this 
is a fair agreement and a good sign of the consistency be- 
tween the present, completely independent estimation of 
Z and those obtained by direct spectroscopic observations 
(of other Cepheid variables). 

Although Fig. 4 shows that the exact DMqx = DM12 
condition is still not satisfied for Z — 0.003, and a some- 
what lower Z would be more appropriate, the differ- 
ence is within the reasonable error limit, and therefore, 
in the estimation of DM, we use the result obtained at 
Z = 0.003. When weighted by the number of stars, we get 
19.05 mag for the average rounded distance modulus. By 
using the lower T e g scale of A99, a similar match is found 
at Z — 0.003 between DMqi and DM\2. In this case we 
get DM = 18.90 mag. Due to the large number of vari- 
ables, the formal statistical errors are very small in both 
cases: cdm = 0.017 mag. Much more significant sources of 
errors are the various systematic effects and zero point am- 
biguities. Table 1 summarizes the changes in the distance 
moduli caused by the most significant potential sources of 
these kinds of errors. It is difficult to assess the size of the 
systematic errors in the various quantities. The numbers 
entering in the table are our best guesses on the 3u errors 
(for notation simplicity we used positive changes every- 
where). The assumed ambiguity in T e g is based on the 
difference between A99 and BLG94. We think that this 
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Fig. 4. As in Fig. 3, but for different heavy element abun- 
dances. 

Table 1. Systematic errors in the distance moduli at 
Z = 0.003 



Quantity 




ADMoi 


ADM 12 


AlogTrff 


= +0.01 


+0.22 


+0.11 


AU 


= +0.02 


-0.06 


-0.03 


AI C 


= +0.02 


+0.10 


+0.05 


AEb-v 


= +0.05 


+0.15 


+0.02 


AP 


= +0.01 


-0.00 




APi/Po 


= +0.002 


-0.15 




APi 


= +0.01 




-0.04 


APa/Px 


= +0.002 




-0.01 



Note: The period errors refer to the model values. 

is a generous overestimation of the true error, because of 
the ambiguities mentioned in Sect. 2 and because of the 
smaller difference obtained in a comparison with the other 
current scale of BLG98 (see Sect. 2). 

The following conclusions can be drawn from the table: 

(a) FU/FO variables are (much) more sensitive to the sys- 
tematic changes than FO/SO variables. 

(b) The most serious source of error is the zero point am- 
biguity in the T e g scale. 

(c) Considering only the FO/SO variables, which domi- 
nate the value of the average distance modulus, and 
assuming that the various errors are independent, we 
get ±0.13 mag for the total estimated systematic error. 

4. Conclusions 

By using the relative distance of 0.51 mag determined 
by Udalski et al. (1999b), the present determination of 
the SMC distance leads to an LMC distance modulus 



of 18.54 mag. This value is magically close to the value 
of 18.53 mag, derived from the application of the same 
method to Galactic double-mode RR Lyrae stars and us- 
ing the relative distance of a few LMC globular clusters 
(Kovacs 2000). Considering that this result was derived on 
a completely different data set, we think that the agree- 
ment is remarkable. This result suggest that the present 
models and input physics are more compatible with the 
observations than implied by Buchler et al. (1996) from 
their beat and bump Cepheid studies. 

As given in Table 1, the most important source of am- 
biguity in this method is the potential error in the zero 
point of the temperature scale. Even if we consider this 
ambiguity, it is not possible to lower the distance mod- 
ulus by more than « 0.13 mag. This emphasizes further 
the contradiction between this 'long' and other 'short' dis- 
tances, obtained e.g., by statistical parallax and red clump 
methods (Udalski et al. 1999b, see however Romaniello et 
al. 2000). 
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